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ABSTRACT 

At low redshift, there are fundamental correlations between the mass of supermassive black holes 
(Mbh) and the mass (Albuige) and luminosity of the host galaxy bulge. We investigate the same 
relation at z > 1. Using virial mass estimates for 11 quasars at z > 2 to measure their black hole mass, 
we find that black holes at high-z fall nearly on the same A^bh versus i?-band magnitude (-Mr) relation 
(to ~ 0.3 mag) as low-redshift active and inactive galaxies, without making any correction for luminosity 
evolution. Using a set of conservative assumptions about the host galaxy stellar population, we show 
that at z > 2 (10 Gyr ago) the ratio of A^bh/AW bulge was 3-6 times larger than today. Barring unknown 
systematic errors on the measurement of A^bh, we also rule out scenarios in which moderately luminous 
quasar hosts at z > 2 were fully formed bulges that passively faded to the present epoch. On the other 
hand, 5 quasar hosts at z « 1 are consistent with current day Mbh-Mr relationship after taking into 
account evolution, appropriate for that of E/S0 galaxies. Therefore, z k, 1 host galaxies appear to fit the 
hypothesis they are fully formed early-type galaxies. We also find that most quasar hosts with absolute 
magnitudes brighter than Mr = —23 cannot fade below L» galaxies today, regardless of their stellar 
population makeup, because their black hole masses are too high and they must arrive at the local 
■Mbh-Mr relationship by z = 0. 

Subject headings: galaxies: evolution — galaxies: quasars — galaxies: fundamental parameters — 
galaxies: structure — galaxies: bulges 



1. INTRODUCTION 

There are currently ~ 40 galaxies in which a supermas- 
sive black hole (BH) has been directly measured using stel- 
lar and gas dynamics (e.g., reviews by Barth 2004 and Ko- 
rmendy 2004). With a large number of secure detections, 
several correlations have been discovered between the mass 
of the BHs and the global properties of the galaxy bulge, 
such as its stellar velocity dispersion (ct*), mass ( Al bulge) > 
and luminosity (e.g., Kormendy & Richstone 1995; Magor- 
rian et al. 1998; Ho 1999; Gebhardt et al. 2000a; Ferrarese 
& Merritt 2000; Kormendy & Gebhardt 2001, Haring & 
Rix 2004). What makes the correlations even more strik- 
ing is the relatively small amount of intrinsic scatter; the 
relation between A^bh and c* (Gebhardt et al. 2000a, and 
Ferrarese & Merritt 2000) has only a scatter of < 2.2 (0.3 
dex; Tremaine et al. 2002) in A4bh- In the original dis- 
covery, the relationship between A^bh and the bulge lumi- 
nosity appeared weaker, with a scatter of 0.5 dex (Magor- 
rian et al. 1998, Kormendy & Gebhardt 2001). When 
the bulge luminosity (ibuige) is translated into mass, one 
finds that the ratio of the BH mass to the bulge mass 
is A4 B H/A4bui g c « 0.0012 (Kormendyfe Gebhardt 2001, 
Merritt & Ferrarese 2001, and McLure & Dunlop 2002). 



These remarkable findings leave little room to doubt that 
the formation and evolution of BHs and galaxy bulges are 
closely tied. Understanding why provides fundamental in- 
sights toward a coherent understanding of galaxy forma- 
tion and evolution in general. 

To study the coevolution of black holes with galaxy 
bulges, there is a great deal of interest in investigating the 
correlations at higher redshifts. In this study, we explore 
the A^bh versus the i?-band luminosity (Mr) relation of 
the galaxy bulge out to z w 2. The Mbk-Mr correlation 
has been studied by McLure & Dunlop (2002), and Bettoni 
et al. (2003) for z < 1 galaxies, and is an extension of the 
studies pioneered by Kormendy & Richstone (1995) and 
Magorrian et al. (1998). Since Kormendy & Richstone's 
(1995) original finding, which was based on £?-band images 
of nearby galaxies, there have been many efforts to sharpen 
the correlation using larger samples and/or redder pass- 
bands (Magorrian et al. 1998; Ho 1999; Kormendy & Geb- 
hardt 2001; Laor 2001; Merritt & Ferrarese 2001; McLure 
& Dunlop 2002; Bettoni et al. 2003; Marconi k, Hunt 2003, 
Ivanov & Alonso-Herrero 2003). The RMS scatter of the 
correlation, when both early- and late-type galaxies with 
A^bh measurements are taken into account, is fairly large 
(~ 0.5 dex; Kormendy & Gebhardt 2001; Marconi & Hunt 
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2003) . When attention is focused on just the early-type 
galaxies, several studies find the scatter in the A^BH-bulge 
luminosity relation drops to 0.3-0.45 dex (Erwin, Graham, 
& Caon 2002, McLure & Dunlop 2002; Bettoni et al. 2003; 
Marconi & Hunt 2003). While it is not clear what causes 
the larger scatter in late-type galaxies, one of the compli- 
cations is that it is not always straightforward to obtain 
bulge luminosities, which is hard to uniquely separate from 
the disk in late-type galaxies when the bulge does not have 
a classical de Vaucouleurs profile (e.g. Carollo 1999). 

While measuring the bulge luminosity of normal (i.e. 
inactive) galaxies is feasible out to high z, measuring the 
mass of their central black holes accurately is considerably 
more difficult. Currently, the most robust method for mea- 
suring A^bh in normal galaxies is through modeling the 
stellar and gas dynamics (e.g., Barth 2004 and Kormendy 

2004) . Doing so requires an exquisite spatial resolution on 
the order of the BH sphere of influence, r sp h = GAIbh/c*, 
the radius at which the Keplerian orbital velocity of the 
stars due to the influence of the central BH is comparable 
to the stellar velocity dispersion of the bulge. This is fea- 
sible only for the nearest, and often relatively luminous, 
galaxies. The A^bh-c* relation in nearby galaxies offers an 
easier way to infer A^bh- However, despite the potential 
for extending this technique widely to infer Mbh to study 
bulge and BH evolution, in practice the technique is time 
consuming and has practical challenges: cr» is difficult to 
measure in a variety of circumstances, including low-mass 
and low-surface brightness galaxies, distant galaxies, and 
galaxies in which the active nucleus substantially dilutes 
the starlight (Greene & Ho 2005a). More fundamentally, 
to study the co-evolution of BHs with galaxy bulges, it is 
desirable to measure A^bh using a technique that is more 
directly associated with the influence of a black hole, es- 
pecially at high z, rather than one which is tied to the 
bulge. 

On the other hand, active galactic nuclei (AGNs), whose 
energetic output is a direct manifestation of supermas- 
sive black holes, are natural and promising candidates for 
measuring A^bh, using a technique known as the "virial" 
technique (e.g. see Ho 1999; Wandel, Peterson, & Malkan 
1999). This technique is based on measuring AGN broad- 
line region (BLR) sizes and linewidths that have been 
subjected to reverberation mapping (Blandford & McKee 
1982; Peterson 1993). These simple A^bh estimates ap- 
pear to obey the A^bh-c* relation of inactive galaxies at 
low redshift (z < 0.05) (Gebhardt et al. 2000b; Ferrarese 
et al. 2001; Nelson et al. 2004, Onken et al. 2004; Barth, 
Greene, & Ho 2005; Green & Ho 2005b), which provides 
an important cross-check on the validity of the AGN virial 
masses. For a sample of 34 low- z Seyfert 1 nuclei and 
quasars with reverberation mapping data for the H/3 line, 
Kaspi et al. (2000) showed that the BLR size correlates 
strongly with the AGN continuum luminosity. This pro- 
vides a convenient shortcut to estimate BH masses essen- 
tially for any (unbeamed) broad-line AGN with a rest- 
frame optical spectrum. It offers a potential to extend the 
measurement of A^bh out to a large sample of objects, 
spanning a wide range in redshift. Since then, this virial 
mass measurement technique has been bootstrapped to 
higher redshifts using the ultraviolet lines C iv A1549 A 
(Vestergaard 2002) and Mg n A2798 A (McLure & Jarvis 



2002) . The calibration sample, however, is small, currently 
based on < 60 objects. 

Using the virial mass technique, the measurement of su- 
permassive BHs has now been extended out to large red- 
shifts in a large number of quasars. Their numerical abun- 
dance at high redshift (z > 1) and their high luminosities 
make them ideal primary tracers for BH evolution, thus 
their host galaxies prime candidates for studying galaxy 
evolution. The trade-off for having an easy proxy for mea- 
suring A^bh in quasars is that detecting quasar host galax- 
ies is a much harder task compared with studying inactive 
galaxies at high z (i.e., z > 1). The difficulty in extract- 
ing robust parameters for the host galaxy in the glare of 
a luminous quasar is well known. It is worth pointing out 
that even with the exquisite 0'.'05 resolution images of the 
Hubble Space Telescope (HST) in the optical, the task re- 
mains nontrivial because of the small host galaxy size and 
its low surface brightness, which decreases as (1 + z) A . De- 
spite technical challenges, careful observations are detect- 
ing more and more quasar hosts at high z, even from the 
ground, using deep imaging and adaptive optics (Fynbo, 
Burud, & M0ller 2002; Lacy et al. 2002; Hutchings 2003; 
Sanchez & Gonzalez-Serrano 2003, Kuhlbrodt et al. 2005). 

In this study we explore the Mbh vs. Mr relation- 
ship for quasar host galaxies out to z > 1. We use the 
virial technique to estimate A4bh in quasars, combined 
with i?-band luminosity of the bulge of quasar host galax- 
ies, inferred from HST studies of Kukula et al. (2001), 
Ridgway et al. (2001) and Peng et al. (2005). At low z 
quasars brighter than My = —23.5 primarily live in el- 
liptical galaxies (McLeod & Rieke 1995 and Dunlop et al. 

2003) , and there are hints that indicate hosts at z > 1 
may be early-types as well. In particular, Kukula et al. 
(2001) and Sanchez & Gonzalez-Serrano (2003) find that 
the luminosities and scale lengths of galaxies are consis- 
tent with a passively evolving scenario, and, moreover, 
that they follow the Kormendy relation (the correlation 
between half-light surface brightness [i X j 2 and effective ra- 
dius ri/ 2 ). Therefore, in this study, we also explore the 
evolution in the Mb^-Mr relationship by testing the null- 
hypothesis that the host galaxies evolve like E/S0 galaxies 
from a formation redshift of Zf = 5. 

We note that at low redshift, the scatter in the AIbh- 
Mr relationship for elliptical galaxies (McLure & Dunlop 
2002 and Bettoni et al. 2003) is small enough that when ex- 
tended to high redshift it can provide a sensitive and useful 
probe of the coevolution of BH masses with bulges. Black 
holes are a useful constraint for understanding galaxy evo- 
lution because A(bh only increases, which also makes it 
a way to age-date a galaxy. Coupled to Lbuige, galaxies 
at any given redshift are restricted in how much they can 
fade in the two-dimensional diagram of Mbh vs. Lbuigc- 
They cannot freely roam the A^BH-^tmigc diagram, lim- 
ited on the one hand by the monotonicity of AIbh, on 
the other by the requirement that they evolve toward the 
tight local A^bh vs. Lbuigc relation at z = 0, and addi- 
tionally by the local galaxy luminosity function and the 
mass function of supermassive black holes. Thus, at a 
given bulge mass, luminosity, and redshift, galaxies may 
follow a range of possible merger and luminosity evolution 
paths that are initially wide but progressively narrow as 
they evolve towards the tight local relations of A^bh ver- 
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sus luminosity. Generalizing the idea into n-dimensional 
space, by measuring the Mbh vs. luminosity relation in 
several passbands, this "phase-diagram" may indicate the 
state of galaxies at high redshift relative to today. This 
may ultimately be a powerful, yet simple way to visual- 
ize and understand the complex paths of galaxy formation 
and evolution. The technique requires calibrating A4bh 
vs. luminosity at many passbands for local galaxies. 

Our discussion below is structured in the following or- 
der. To obtain the A4bh~Mr relation for quasar host 
galaxies, one of the first steps is to make fc-correction to the 
host galaxy photometry (§ 2). In Section 3, we will then 
estimate the BH masses using the virial technique, com- 
pare them to a cruder estimate using the AlBH-himinosity 
relation for quasars determined by Peterson et al. (2004), 
and discuss the Eddington efficiencies of the quasars in the 
sample, using black hole masses estimated from the virial 
technique. Section 4 presents the results, and we conclude 
with brief remarks (§ 5). In this study, we use a standard 
cosmology with Ho = 70 km -1 s _1 Mpc -1 , O m = 0.3, and 
Q A = 0.7. 

2. DATA AND BULGE fc-CORRECTION 

We assemble a set of data on quasar host galaxies pub- 
lished by Kukula et al. (2001), Ridgway et al. (2001), and 
Peng et al. (2005) observed using tfST/NICMOS, summa- 
rized in Table 1. All the magnitudes shown in the Table 
have been corrected for Galactic extinction according to 
Schlegcl, Finkbciner, & Davis (1998), but not internal ex- 
tinction. The Kukula et al. (2001) sample consists of 18 
quasars at z ~ 1 and z « 2, equally divided between radio- 
loud quasars (RLQ) and radio-quiet quasars (RQQ), while 
the Ridgway et al. (2001) sample consists of 5 RQQs at 
1.5<z<2.8. 

The Peng et al. (2005) study has one additional object, 
obtained using the gravitational lensing technique of the 
CASTLES survey (Kochanek et al. 1999). In that study, 
they find that the host galaxy has been stretched out into 
a luminous Einstein ring. After deprojection, where the 
lensing galaxy mass is assumed to be a singular isother- 
mal ellipsoid, they recovered the intrinsic light distribution 
of the host galaxy (e.g. size, luminosity, axis ratio, etc.). 
The host galaxy Einstein ring in CTQ 414 is only mod- 
estly merged with the lensing galaxy light profile. Even 
so, because the shape of the host galaxy Einstein ring 
is so dramatically different azimuthally from the lensing 
galaxy, it is fairly easy to disentangle the two. They did 
so by modeling the lens deflection model, light profile of 
the lens, the host, and the quasar images - all simultane- 
ously. In so doing, they find that the uncertainty of the 
host luminosity is roughly 0.2 mag. As we shall see, the 
result for this object would only start to be important if 
it were one magnitude brighter than our number. Even 
so, individually, no one object will significantly affect the 
overall conclusion. Thus this object alone serves mostly as 
a consistency check, as it is the most well- resolved object. 

The data from all samples were observed in either NIC- 
MOS F110M (~ J band) or F160W/F165M (~ H band) 
filters so that they roughly correspond to rest-frame V, 
where the contrast between the host and the quasar is fa- 
vorable for detecting the host galaxy. Altogether, 24 host 
galaxies are used in this study, which is the entire col- 



lection of data in the literature published specifically to 
study quasar host galaxies detected using NICMOS above 
z > 1. There are other host galaxies found at z > 2 pub- 
lished using WFPC2 and ground-based data (e.g., Aretx- 
aga, Terlevich, & Boyle 1998; Lehnert et al. 1999; Hutch- 
ings et al. 2002: Hutchings 2003). We exclude these be- 
cause the observations were made in the rest-frame UV, 
which is very sensitive to star formation and extinction 
internal to the host galaxy. Therefore, fc-correcting to the 
rest-frame R filter is highly dependent on assumptions of 
the spectral energy distribution (SED) of the galaxy. We 
also exclude ground-based data of RLQ hosts published by 
Falomo, Kotilainen & Treves (2001), Lacy et al. (2002), 
Sanchez & Gonzalez-Serrano (2003), and Kuhlbrodt et al. 
2005, to keep the sample uniform and to avoid potential 
uncertainties with the quasar-host separation relative to 
data obtained uniformly with HST. 

To infer the i?-band absolute magnitude (Mr) of the 
quasar host galaxy bulge from the raw values published 
in Kukula et al. (2001), Ridgway et al. (2001), and 
Peng et al. (2005), we make several corrections. We note 
that the published host galaxy fluxes primarily refer to 
the bulge because the images are too shallow to detect 
the diffuse components of extended disks if they exist, ex- 
cept for CTQ 414 (Peng et al. 2005). For that object, 
we take the bulge luminosity only. We point out that 
since the local black hole versus bulge relation excludes 
the disk component, as we shall see, a proper removal of a 
disk would further strength our conclusions by making the 
bulges fainter. Ridgway et al. (2001) compared their host 
galaxy detection with simulations, from which they deter- 
mined the host galaxy flux as a function of the morpho- 
logical type. Thus, for their sample, to obtain the if-band 
bulge luminosity shown in Table 1, we apply an aperture 
correction (their Table 4) to their l'.'Ol aperture photome- 
try (their Table 2), and a correction that is dependent on 
assumptions about the morphology of the host (their Table 
5), appropriate for radio galaxies, which have the largest 
corrections in their study. Finally, to convert published 
photometry into a rest-frame, standard Cousins, i?-band 
flux, we follow Hogg et al. (2004) in computing detailed 
fc-corrections. We transform the HST Vega-based magni- 
tudes to i?-band by using galaxy templates observed by 
Coleman, Wu, & Weedman (1980), multiplied by the ap- 
propriate filter transmission curves. The integrated fluxes 
are normalized to a spectrum of a-Lyr, "observed" in the 
appropriate bandpass. 

Computing /c-corrections to the host galaxy's rest-frame 
i?-band relies on having an SED. We have no color infor- 
mation for most of the host galaxies so we must make an 
assumption. The filters used by Kukula et al. (2001; J 
and H) and Ridgway et al. (2001; H) conveniently cap- 
ture the hosts in their rest-frame V (except for two objects 
at z s=s 2.7 in Ridgway et al. 2001). Therefore, /c-correction 
to rest-frame R is not very sensitive to assumptions of the 
SED for both the quasar or the host galaxies, especially 
as the V-filter lies redward of the rest-frame 4000 A break 
in galaxies. Specifically, the difference in the inferred Mr 
between an early and late- type SED is only ~ 0.1 mag at 
z = 1 and ~ 0.3 mag at z = 2. For the fc-correction we use 
an SED of current-day E/SOs for the host galaxies because 
it is the most conservative assumption with regard to the 
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conclusions we are testing (i.e. weakening it), despite it 
being less sensible than other potential models. All other 
more realistic SEDs at z ~ 2 are likely to be bluer, thus the 
inferred host luminosities would be systematically fainter 
than the Mr values shown in Tabic 2. 

3. BLACK HOLE MASS ESTIMATE AND EDDINGTON RATIO 

3.1. Virial A4bh Measurement 

For high redshifts, the virial technique relies on mea- 
suring the quasar continuum luminosity and the widths 
(FWHM) of the C iv (A1549 A) or Mg n (A2798 A) emis- 
sion line. At z ~ 1 the Mg n emission line is the most 
convenient feature to use because it falls into optical band- 
passes where the atmosphere is transparent and where 
CCDs are sensitive. Likewise, at z > 2 C iv is the line 
of choice. Compared to the H/3 line, C iv is produced 
much closer to the supermassive black hole (by factor of 
~ 0.5). As such there are concerns that the kinematics of 
the BLR, not well understood, may be dominated by non- 
Keplerian motions such as winds and outflows. Concerns 
about the UV emitting clouds being caught up in a disc 
also meant that the width of the emission lines may de- 
pend strongly on the viewing angle. Despite uncertainties 
in the BLR physical structures, the same virial product 
that is used to estimate black hole masses in AGNs us- 
ing H/3 is found to have analogs in the UV with C iv or 
Mg II lines, over the luminosity range present in our AGN 
sample (Vestergaard et al. 2002, McLure & Jarvis 2002, 
Kaspi et al. 2005). The assumption we will make is that 
the same technique applies to quasars at high redshifts. 

The high redshift calibration for C IV was performed by 
Vestergaard (2002) 9 . In viral estimates of A4bh there is 
a geometric normalization factor, /, that depends on the 
structure and kinematics of the BLR, and which normal- 
izes the AGN A4bh-c* to that of the quiescent galaxies. 
This factor was determined empirically by Peterson et al. 
(2004) and Onken et al. (2004), to be (/) = 5.5 for the 
H/3 BLR (under an assumption that au ne — FWHM/2). In 
the absence of a large sample to determine /, the effective 
normalization used in Kaspi et al. (2000) corresponded 
to (/) = 3, so that the new Mbh estimates are a factor 
of 5.5/3 = 1.8 higher than previously estimated in Kaspi 
et al. (2000). With these caveats, M. Vestergaard (2005, 
private communication) provided the following relation to 
estimate A4bh from C iv (A1549 A), appropriate for our 
adopted cosmology: 



region (McLure & Jarvis 2002). Again, we renormalize 
their empirical relation by a factor of / = 1.8 to obtain: 



M 



BH 



4.5xl0 6 



FWHM(C iv) 



1000 km s" 



AL A (1350 A) 



-, 0.53 



10 44 erg s- 



M 

(1) 

where La (1350 A) is the spectral density of the quasar 
continuum at rest-frame 1350 A. The value of the ex- 
ponent (7 = 0.53) and the coefficient (4.5 xlO 6 ) depend 
on the sample used in the fit and the fitting technique 
(BCES [Akritas & Bcrshady 1996] vs. FITEXY [Press et 
al. 1992]). There may be small systematic differences, but 
they are otherwise formally equivalent to within the errors 
(Kaspi et al. 2005). 

For z « 1 quasars Mg 11 (A2798 A) is a good surrogate 
for estimating AIbh that conveniently falls in the visible 

9 She used a cosmology of Hq = 74 km s _1 Mpc -1 , SI a = 0, and 
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(2) 

To obtain the C iv and Mg 11 emission line FWHM in 
quasars with host detections, we searched the literature 
for published spectra. In total, we compile 16 (out of a to- 
tal of 24) objects that have published spectra or estimated 
FWHM values - 11 at z > 2 and 5 at ~ 1. The Ridgway 
et al. (2001) sample comes from a survey published by 
Zitelli et al. (1992), while the Kukula et al.' (2001) sample 
is a composite of different surveys shown in Table 1. For 
CTQ 414, we obtained an uncalibrated spectrum at the 
Magellan telescope. We note that the quasar spectra pub- 
lished in Boyle et al. (1990) used by Kukula et al. (2001) 
sample have very uncertain flux calibrations, and likewise 
for CTQ 414. However, what matters more to our anal- 
ysis is the FWHM values, given that the absolute fluxes 
are determined by photometry (see below). Moreover, be- 
cause the continuum luminosity enters into Eqs. 1 and 2 
at most as the power of 0.53 (0.47 for Mg 11), even a flux 
uncertainty as large as a factor of 2 translates into an un- 
certainty in A^bh of at most 0.2 dex, which will not affect 
our conclusions. To measure the widths of the emission 
lines that do not have published values, two of the au- 
thors independently measured the FWHM manually from 
compiled spectra, in a double-blind test. The two indepen- 
dent measurements of the FWHM agree well. One other 
potential uncertainty in measuring the FWHM is that the 
broad emission lines may be contaminated by Fe 11 emis- 
sion (Vestergaard & Wilkes 2001). Fe 11 contamination is 
more important for H/3 and Mg 11 lines (Vestergaard & 
Wilkes 2001), and much less for C iv. For the Mg 11 line 
we do not have a way to estimate the Fe 11 contribution, 
thus we do not correct for this effect. This is appropriate 
for the calibration performed by McLure & Jarvis (2002), 
which did not remove the Fe 11 contamination as well. All 
of our z > 2 objects use C iv measurements and only our 
z ?y 1 objects rely on Mg 11 lines. 

To estimate the monochromatic luminosity of the 
quasars at La (1350 A) and La (3000 A), wc first k- 
corrected to rcstframe [/-band, centered on 3745 A, using 
a power law SED that fits to apparent magnitudes in H/ J, 
and B/V filters. The same SED is then used to extrapolate 
to 1350 A or 3000 A as required. To convert from [/-band 
magnitude to flux density we use the spectral irradiance 
value of 4.34 x 10~ 9 erg s" 1 cm~ 2 A -1 (Colina, Bohlin, & 
Castelli 1996). There is a small conversion to take the B- 
band magnitudes from literature, based on photographic 
plates (see Table 1), into the standard Vega-magnitude 
system. The magnitudes of the SGP objects from Boyle 
et al. (1990) were based on plate-scanned Bj system ("6- 
band"), zeropointed to the Johnson B magnitude system. 
The conversion to the standard B magnitude has a small 
offset for typical quasar colors: b = B + 0.1 ± 0.05. In 
comparison, the random errors are generally larger, about 
0.2-0.3 mag. Similarly, the magnitude of the MZZ objects 
from Zitelli ct al. (1992) is based on photographic Bj filter 
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system. As the systematic corrections are small and will 
not affect our conclusions, we do not apply a correction. 
We correct the B-band magnitudes for Galactic extinction, 
according to Schlcgcl, Finkbeincr, & Davis (1998). 

This technique for measuring the quasar flux density 
near 1350 A and 3000 A is fairly robust, as the B/V and 
J I H filters cover a long baseline for estimating a powerlaw 
SED. Since the observed B band also falls in the restframe 
UV, there is only a small amount of SED-dependent ex- 
trapolation. We compare A^bh estimates using our SED 
extrapolation with another crude method of directly con- 
verting observed i?-band and U-band fluxes into Mbh- 
Both methods give nearly the same results. Contamina- 
tion by emission lines in the broad-band filter is often neg- 
ligible even in quasars because the bandwidths of the V 
(FWHM = 890 A) and B (FWHM=1000 A) filters are 
significantly larger than the widths of the lines (FWHM 
< 100 A). The impact of such errors on AIbh is negligi- 
ble (< 5%) because the monochromatic luminosity enters 
into Eqs. 1 and 2 as the powers of 0.53 and 0.47, respec- 
tively. All the transformed values are shown in Table 2. 
For CTQ 414 and 4C 45.51 we take the average of A^bh 
estimates from the C iv and Mg n techniques. 

Lastly, we point out that with only two objects, 4C 45.51 
and CTQ 414, that have both C iv and Mg n simultane- 
ously observed it is not possible to check on systematic 
differences in the A^bh estimates: in 4C 45.51, the in- 
ferred black hole mass is comparable between using C iv 
(0.65 x 1O 9 A4 ) and Mg n (0.58 x 10 9 Af©). In the other, 
CTQ 414, C iv yields a factor of 2 larger BH mass than 
Mg II. However, in a separate study using gravitationally 
lensed quasars at z > 1 (Peng et al. in prep), including 4 
other objects that have simultaneous C iv and Mg n obser- 
vations, the two methods appear not to have a systematic 
offset. 

3.2. A4bh- Luminosity (Lq$o) Relation of Quasars 

Using virial estimates of A4bh, it has been pointed out 
that, to first order, BH masses correlate with the unre- 
solved AGN optical/UV luminosity (e.g. see Kaspi et 
al. 2000; Oshlack, Webster, & Whiting 2002; Peterson 
et al. 2004). This correlation is an empirical statement 
that while the distribution of Eddington ratio e in AGNs 
is large, there might be a preferential efficiency over an 
average lifetime once they are "turned on." This typi- 
cal e might also depend on the luminosity range of AGN 
sample under consideration. For nearby AGNs the mass- 
luminosity relation of quasars is shown in Peterson et al. 
(2004) to be 10 : 

log(A4 BH /10 8 M Q ) = -0.12(±0.17)+ 

0.79(±0.09) log(AL A (5100 A)/10 44 ergs" 1 ). (3) 

For objects which we do not have spectral information to 
determine accurate virial BH masses, we can obtain a more 
crude estimate through this AIbh-Lqso relation. We also 
point out that Eq. 3 is based on measurements of nearby 
AGNs, and currently there is little reason to believe that 
the same normalization would apply to luminous quasars 
at z — 2 where the distribution of luminosities, Edding- 
ton ratio, and AGN lifetimes, may be somewhat different 

10 All logarithms are base 10. 



(e.g. see Kollmeier et a. 2005). The goal of this crude 
method is thus only to gauge how a bigger scatter, uncer- 
tainties in the overall BH mass normalization and small 
number statistics, can affect interpretation. And wc em- 
phasize that our final conclusions are based only on virial 
A4bh estimates alone. 

3.3. Eddington Ratios of Quasars 

From A^bh estimates using the virial technique 
and assuming that the bolomctric luminosity Lboi ~ 
9ALa(5100 A), as also assumed by Peterson et al. (2004) 
and Kaspi et al. (2000), we can estimate the Eddington 
ratio, e, for 11 quasars at z « 2, and 4 at z « 1, which have 
published spectra as a check for any unusual objects (e.g. 
e > 1). For z « 1 quasars which lack the required U-band 
magnitude, we assume B — V = 0.3, which corresponds 
to a continuum of f v cx i/~ 05 (Vanden Berk et al. 2001). 
Table 2, Col. 11, shows the Eddington ratio computed 
for 15 quasars that have virial A4bh estimates. All the 
quasars in our sample have Lboi < -^Edd, with the only 
object, 4C 45.51, which is radio loud, near the Edding- 
ton limit. The Eddington ratio spans a wide range from 
0.02 < e < 0.7, with the top two highest efficiencies being 
RLQs and the bottom two being RQQs. This is probably 
due to a selection effect since Kukula et al. (2001) selected 
RLQ and RQQ quasars that are well matched in optical 
luminosities, but required RLQs to be extremely radio lu- 
minous (> 1 Jy sources, or ~ 3 x 10 28 to ~ 10 29 W Hz -1 ). 
Because our sample size is small, and because a number of 
z « 2 RLQs and z w 1 RQQs are missing virial A4bh 
measurements, the issue of a correlation between radio 
loudness and BH mass and e, at a given rcdshift, can- 
not be addressed here. We only note that such claims of 
a correlation are controversial even at low-z: while several 
studies (e.g. Pagani, Falomo, & Treves 2003; Dunlop et al. 
2003; and McLure & Dunlop 2001) find that RLQs tend 
to have more massive A4bh and more powerful quasars 
than RQQs at low z, several other studies find little to 
no correlation (e.g. Ho 2002; Woo & Urry 2002; Oshlack, 
Webster, & Whiting 2002; Gu, Cao, & Jiang 2001). 

3.4. Measurement Uncertainties 

In this study, the robustness of the conclusions relies on 
understanding the sources for both systematic and random 
errors. In light of systematic errors, the best way to min- 
imize potential doubts is to apply all known systematic 
errors (by offsetting the data) towards weakening the final 
conclusions. Specifically, this can be done by using con- 
versions that yield systematically small AIbh an d bright 
host luminosities. We therefore take measures to do so, as 
described below. In this section we summarize the known 
sources of random and systematic errors, some of which 
are collected from previous discussions. 

There are several important sources of random errors 
when trying to interpret the scatter in the plots. In par- 
ticular, the literature measurement of the emission-line 
widths are probably uncertain by ~ 20 — 30% due to man- 
ual measurements and complications from absorption-line 
features and probable Fe n contamination. This produces 
roughly a random scatter of <~ 0.2 — 0.3 dex in AIbh, while 
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lines which are absorbed in the wings may produce un- 
derestimates. The uncertainty in the Mbh measurement 
using the virial technique has about 0.3-0.4 dcx in scatter 
(factor of 2.5; Kaspi et al. 2000, Vestergaard 2002, and 
McLure & Jarvis 2002, Peterson et al. 2004), which is also 
random. The uncertainties of the host galaxy magnitudes 
are better quantified: typically <~ 0.2 — 0.75 mag (Ridg- 
way et al. 2001 and Kukula et al. 2001). /c-correction 
introduces both a systematic bias (by assuming the red- 
dest possible SED) and a random component (the intrinsic 
SED scatter of the hosts). In the Mbh-Lqso method of 
estimating A^bh, one uncertainty is in the intrinsic scat- 
ter in converting optical luminosity, using ALa(5100 A), 
into bolometric luminosity. This is effectively factored into 
the estimate of A4bh from ALa(5100 A), shown in Eq. 3. 
As will be evident, photometric errors (quasar and host 
galaxy luminosities) will not obfuscate the conclusions, be- 
cause they are small compared to the intrinsic scatter in 
empirical relationships between Mbh-Mr, and the virial 
technique in estimating Mbh- And the results are robust 
against even these random sources of error. 

On the other hand, the robustness of the conclusion will 
depend entirely on the treatment of systematic issues. We 
address this issue by applying known errors to the data in a 
way to yield a smaller A^bh and more luminous host galax- 
ies, thereby weakening the conclusions in the outcome. For 
instance, in computing the host galaxy Mr, the main sys- 
tematic error comes from the quasar removal. In Ridgway 
et al. (2001) they find after running simulations that the 
amount of over-subtraction of the host galaxy flux depends 
on the host morphology. Thus, we choose a morphology 
correction that results in our hosts being typically 0.3-0.5 
mag brighter than quoted in their study. In Kukula et 
al. (2001), they perform 2-dimensional image decomposi- 
tion, using a de Vaucouleurs model profile to extrapolate 
the host both into the core and out to the wings. Fortu- 
nately, assuming a de Vaucouleurs profile often tends to 
overestimate host fluxes if the intrinsic profile is shallower, 
e.g. an exponential. Furthermore, they hold fixed the ef- 
fective radii to 5-10 kpc in the fit - large compared with 
even typical quasar hosts at z > 2 (i.e. measured through 
gravitational lensing, Peng et al. in prep.), and Lyman 
break galaxies (Trujillo et al. 2005). We thus feel confi- 
dent about taking their host measurements at face value 
to be reliable upper limits. In converting to restframe R- 
magnitude, there are uncertainties also in the fc-correction: 
we again decide to err on the side of overestimating the 
host luminosity by assuming a red SED, typical of present- 
day E/S0 galaxies. The difference between a red SED and 
a typical Lyman-break galaxy SED is only about 0.3 mag- 
nitudes in restframe R out at z = 2, when observed in the 
iJ-band. In addition, most of the NICMOS images are not 
sufficiently deep to ascertain whether the hosts are bulge 
or disk dominated; we attribute all the host galaxy light 
to the bulge, which is also a conservative assumption. 

The systematic errors for estimating AIbh are managed 
in the following way. In particular, there is some uncer- 
tainty in the luminosity dependence (expressed via the 
exponent, 7, in Eqs. 1 and 2) of the broad line radius 
in using the virial mass estimators. Depending on the 
line fitting technique used, and which sub-sample is being 
studied (e.g. high-luminosity vs. low-luminosity) 7 ranges 



from 0.5 < 7 < 0.74, although a lower 7 is preferred in the 
newest calibrations for the UV broad line regions (Kaspi 
et al. 2005). For the C iv virial A4bh estimate, the new 
calibration uses a value of 7 = 0.53 (M. Vestergaard, pri- 
vate communication), while 7 = 0.47 for Mg 11 (McLure 
& Jarvis 2002). Also, in the virial A^bh estimate (§ 3.1) 
there is a potential systematic error in the "/-factor" that 
normalizes the virial A^bh to the A4bh-c*- This normal- 
ization of / = 5.5, which increases A4bh by a factor of 
1.8 compared to Kaspi et al. (2000), is based on a mod- 
est sample size of 16 objects (Onken et al. 2004). Since 
this is the best determined value currently we will utilize 
it. However, leaving it out weakens the result by 0.5 mag- 
nitude in Mr, which, by coincidence, corresponds closely 
to results using the mass-luminosity relation to estimate 
A^bh- Even this measure would leave the qualitative con- 
clusions unchanged, especially in light of other steps we 
have taken to weaken the conclusions; we will present both 
techniques side-by-side to emphasize this fact. In estimat- 
ing BH masses using the Mbh-Lqso technique, Peterson 
et al. (2004) discuss in more depth the associated uncer- 
tainties. We note that the Y-intercept, —0.12, of Eq. 3 
has a significant uncertainty of ±0.17, potentially causing 
a small normalization difference in mass estimates between 
this and the virial technique. We also do not expect the 
normalization for low redshift AGNs be identical to high 
rcdshift objects because of different selection, and perhaps 
physical, issues. 

There is one particular object, the RLQ 4C 45.51, which 
deviates significantly from others and it is worthwhile to 
discuss the relevant uncertainties. The uncertainty in the 
BH mass is not known, because the uncertainty in measur- 
ing the FWHM is not known. If one assumes the FWHM 
is in error by 50%, an unlikely high error, the error in esti- 
mating A^bh is ~ 0.4 dex. The potential that the quasar 
is beamed adds an unknown systematic error into the es- 
timate of A^bh- For the host galaxy, Kukula et al. (2001) 
quotes a generic errorbar for the ensemble of 0.75 mag, 
barring small uncertainties in /c-correction. All the known 
errors added constructively cannot explain its large depar- 
ture from the other data points. Other indicators, such as 
the high Eddington ratio (e ~ 0.7) and the extreme radio 
loudness (~ 3 x 10 28 to - 10 29 W Hz" 1 from 151 MHz to 
37 GHz, e.g. When et al. 1992; Hales, Baldwin & Warner 
1993) point to the possibility that 4C 45.51 may be an 
exceptional object. Thus we will exclude this object from 
our discussions and analysis below. 

Lastly, our conclusions rest crucially on the applicabil- 
ity of local virial calibrations to objects out at consider- 
ably higher redshift, z = 2, and that they apply for RLQs 
whose emission may be beamed towards us (Oshlack, Web- 
ster, & Whiting 2002; Greene & Ho 2005c). The propriety 
of applying these calibrations to luminous quasars is more 
rigorously discussed, e.g. in Netzer (2003) and Vestergaard 
(2004) ; we point out that the luminosity range of quasars 
in this sample is spanned by those used to calibrate the 
virial relationships locally. 

In summary, to mitigate against systematic errors that 
may complicate the interpretation, we have added known 
systematic offset corrections to the data that would most 
weaken the conclusions. Therefore, any remaining devia- 
tion that we subsequently observe might be a lower bound 
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on the true difference if all the sources of systematic errors 
are identified as discussed here. 

4. RESULTS 

4.1. Mbh-Mk Relationship for z > 1 Quasar Hosts 

Figure 1 shows the Mbh-Mr relationship derived for 
high-z quasars (open symbols) based on the virial mass 
measurements (Fig. la) or on the mass versus quasar lu- 
minosity relation (Fig. 1 6) . The round solid points ( u z « 
norm, gal.") are local calibrations of normal elliptical 
galaxies obtained in various studies. In particular, the 
Mbh measurements for z « sample are the revised val- 
ues given in Trcmaine et al. (2002), while their photome- 
try is compiled from different sources and reprocessed into 
i?-band values by Bettoni et al. (2003). There are many 
measurements of other galaxy types that can be (and have 
been) placed onto this plot, such as radio galaxies (Bet- 
toni et al. 2003) and low- z quasar and Seyfert hosts (e.g., 
McLure & Dunlop 2002). All of them essentially fall along 
the solid line shown, and have comparable scatter (McLure 
& Dunlop 2002; Bettoni et al. 2003). The solid lines shown 
in Figures 1 a and 1 b are based on a fit to the normal galax- 
ies data points (round solid) performed by Bettoni et al. 
(2003). Converting their cosmology (Hq = 50 km s _1 
Mpc -1 ) into ours, the line has the following relation: 

log (Mbh/Mq) = -0.50(±0.06)M fl - 2.70(±1.35). (4) 

Consistency in the Mbh-Mr Relation Between z « and 
z > 2. Figures la and 16 show a comparison of the tech- 
niques used to estimate BH masses out at high z, illustrat- 
ing curiously that most of the z > 2 host galaxies already 
lie near the Mbh-Mr relationship at z w 0. The z~l 
quasar hosts appear to lie noticeably below the relation. 
However, as we have not yet accounted for luminosity fad- 
ing, it is not yet obvious what is to be expected. Most of 
the attention should be focused on Figure la because the 
estimates of A4bh are more robust. Figure 16 merely con- 
firms the fact that even with a larger scatter, and possibly 
a systematic normalization offset in BH mass, the result 
does not change significantly. Although the points using 
the crude estimate agree qualitatively with the virial es- 
timates there are some notable exceptions that bring the 
average in Figure lb considerably lower. The most dis- 
crepant points are the two extreme radio sources objects 
4C 45.51 and B2 2156+29 from Kukula et al. (2001), at 
least one of which has high Eddington ratio. In Figure la, 
the correlation seen for z > 1 objects is probably real and 
not due to strong selection effects in detecting the host 
galaxy, despite well-known challenges in extracting host 
galaxies from underneath luminous quasars. While one 
does expect selection bias to hamper finding hosts lying 
to the left of the sample points, this appears not to be 
a serious problem because 23 host are detected out of 24 
targets in Kukula et al. (2001) and Ridgway et al. (2001). 
On the other hand, if host galaxies exist to the right of the 
sample points, they would be easier to detect. Thus, the 
fact that there is not a larger and a more uniform scatter 
to the right suggests that the Mbh-Mh relation has al- 
ready been established as early as z — 2, at least in RQQ 
hosts. 



To further interpret Figures la and 16, it is useful to 
keep in mind that regardless of where the high-z quasar 
host galaxies lie, by z — they must scatter around the 
solid line representing the relationship for z = galaxies. 
Thus, it is surprising that ignoring any luminosity evo- 
lution, the Mbh-Mr relation of quasar host galaxies at 
z > 2 already falls near the same relation as local galax- 
ies. To quantify the similarity of z > 2 AGN hosts to the 
local relation we fit a line to the RQQ objects (shown as 
dotted lines in Fig. 1), holding the slope fixed to the value 
of —0.5 determined for low- z galaxies (Equation 4). As 
a caution we leave RLQ objects out of the fit because of 
possible issues with beaming and because of their higher 
Eddington ratios, as discussed in Section 3.3, even though 
most objects do not deviate far from RQQs. There is only 
one object which is radio loud in Figure la (4C 45.51) and 
4 objects in Figure lb- two of which have the highest host 
luminosity in the sample, and are the farthest outliers. 

In Figure la, the horizontal offset between the dotted 
and solid lines is —0.3 mag when RLQ, 4C 45.51 is ex- 
cluded. If included, 4C 45.51 would single-handedly shift 
the dotted line further to the right by 0.3 magnitudes, even 
though it is only one out of 11 objects. The conclusions are 
largely immune to this shift, as evidenced later when evo- 
lution models are factored in. The average shift in Figure 
la is also somewhat biased by a single RQQ, MZZ 9592 
(green triangle near Mr« —25). This host is the second 
highest rcdshift object (z = 2.71) in the sample, observed 
in restframe i3-band. Compared to a blue powerlaw SED 
{v~ x ) of coeval Lyman-break galaxies, this object has k- 
correction value brighter by 1 magnitude (as opposed to 
~ 0.3 mag for other hosts). In light of this and our use of 
the E/S0 SED for fc-correction, if the hosts are as blue as 
Lyman break galaxies, the agreement between z w 2 and 
the local Mbh-Mr relation would be even more striking, 
bringing the two closer by at least 0.3 magnitudes, bridg- 
ing any difference between high z bulges and z = normal 
galaxies. 

In Figure 16, the offset in the dashed line is larger, —0.6 
mag, which indicates a normalization difference between 
the two estimates of A^bh- This line also excludes 4 out 
of 10 objects, which are RLQ hosts; including them would 
shift the dashed line to the right by 0.2 magnitudes. 

A-^BH/Mbuige Ratio at z w 2. If the local Mbb~Mr re- 
lation more fundamentally reflects a relationship between 
A^bh and A4 bulge, rather than with the bulge luminosity, 
then the fact that z > 2 hosts lie near the same relationship 
as z rts data is unexpected even before applying a correc- 
tion for luminosity evolution. About half of z « 2 hosts 
are nearly on the z ps relation - more, if the SEDs of 
the hosts are bluer (hence fainter) than we have assumed. 
If the hosts starved their black holes (i.e. constant Mbh) 
and there has been no galaxy merger since z ~ 2, this 
implies that the bulge mass was lower at high z: as the 
host faded, more stars would have to be formed in concert 
to build up the bulge mass (luminosity), allowing the host 
galaxy to remain roughly fixed in the A^BH-^i? diagram. 
We can constrain the A^BH/A^buigo ratio at z « 2 if we as- 
sume an evolutionary scenario that the RQQ host galaxies 
at z > 2 are fully formed ellipticals. The offset between 
the dotted and solid lines (Figure la) allows a fading of 0.3 
mag on average since z > 2 (less, if A^bh also increases or 
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if the bulge SED is blue). Such a small amount of dimming 
is perhaps surprising, and we will discuss several possible 
evolutionary scenarios below. We show later that under 
the premise that the hosts evolve like E/SO galaxies, simple 
evolution models predict fading by factors 1-2 magnitudes 
between z = 2 and z = (Figure 2). In this scenario, Fig- 
ure 1 then indicates that the AIbh/ Af bulge ratio at z = 2 
is higher than locally by a factor of 3-6. Rix et al. (1999, 
2001), through studying the hosts of gravitationally lensed 
quasars, first noted that unless all the quasars are radiat- 
ing near (e > 0.5) the Eddington limit, the A^BH/A^buige 
ratio must be higher at high-z than today. Ridgway et al. 
(2001) also came to this same conclusion using the same 
assumption. Here, the argument is made more stringent 
because the black hole mass is measured simultaneously 
with the host galaxy light, and the fact that Mbh can 
only grow with age. 

Deviations from the Mbh-Mr Relation. The largest 
departure from the dotted line in Figure la is 4C 45.51 at 
z = 2, and the two largest in Figure lb are also radio loud 
quasar hosts, one of which is 4C 45.51. Taken at face value, 
its deviation is so large that the apparent correlation seen 
for most z « 2 hosts in Figure la may not be clear-cut 
for all objects, if beaming is not an important factor to 
consider when estimating A^bh- However, it is important 
to keep in mind that this may have nothing to do with the 
issue of RLQ vs. RQQ hosts specifically. Rather, it could 
be a statement about massive vs. less massive galaxies in 
general. The object 4C 45.51 is a radio source whose ex- 
treme radio emissions would likely require an unusually 
massive BH to produce. If the A^BH-Afbuige relation con- 
tinues to hold at high-z, then a selection on extreme AIbh 
(due to selecting on powerful radio sources) would then 
more efficiently pick out massive hosts, even if RLQ and 
RQQs are hosted by the same underlying parent distribu- 
tion of galaxies. The location of 4C 45.51 in Figure la, 
especially, may show that massive, thus rare galaxies may 
be closer to their end-stage of evolution at z s=a 2 than less 
massive ones. 

On the other hand extreme radio loudness can itself 
indicate beaming, which might cause A^bh to be under- 
estimated in RLQs. If the RLQ emission is in fact beamed, 
the continuum luminosity alone would lead to an overes- 
timate of AIbh by a factor proportional to <~ 7l, where 
"fL now refers to the Lorentz factor, and a proportionality 
constant that depends on the viewing angle. For 4C 45.51, 
it may be that 7l > 1 because of its near-Eddington lumi- 
nosity. On the other hand, if the kinematics of the BLR are 
disc-like and the axis of the jet is aligned with the axis of 
the disc, we expect the FWHM of the broad emission lines 
in beamed objects to appear systematically narrower than 
their intrinsic width (Wills & Browne 1986). Since A4bh 
scales as FWHM 2 , the inclination effect may counteract 
the Doppler boosting. Thus, the net effect on measur- 
ing AfsH depends on better knowing the BLR inclination 
and the intrinsic 7l- For RLQs with low inclination and 
low Doppler boosting, we probably underestimate AJbh- 
However, it is difficult to draw a strong conclusion without 
a larger sample of objects. 

Finally, as we already pointed out but worth emphasiz- 
ing, the horizontal shift in Mr of the points in Figure 1 
depends rather little on the E/S0 SED we assumed, and 



uncertainty about the potential of a disk component in the 
host. If we had assumed either a bluer SED, or assume 
that some fraction of the bulge light belongs to the disk, 
they would both necessarily shift the points to the left, 
further strengthening the conclusion that hosts at z = 2 
already lie close to the local Mbh-Mr relation. 

4.2. Host Galaxy Luminosity Evolution 

The near coincidence in the Mbr-Mr relation between 
z > 2 RQQ hosts and the present epoch must reflect sev- 
eral compensating effects, as time has not stood still since 
z ~ 2. Therefore, it places interesting constraints how the 
host galaxies can evolve since the Universe was 2-3 Gyr 
old. At low z (z < 1), McLeod & Rieke (1995) and Dun- 
lop et al. (2003) show that most very luminous quasars 
(My < —23.5 mag) live in ellipticals, with the variety of 
host galaxies increasing towards lower-luminosity AGNs. 
At intermediate redshifts (1 < z < 1.5), several studies 
(Kukula et al. 2001, Sanchez & Gonzalez-Serrano 2003) 
also find that the hosts may be preferentially early-type. 
At z > 1.5, host galaxies are poorly resolved, but there is 
circumstantial evidence to suggest that they are also pri- 
marily ellipticals, fully formed. Hence, we test the quite 
natural null-hypothesis that the quasar host galaxies had 
fully formed stellar bodies and have merely faded since 
z ~ 2. Alternatively, the hypothesis can be phrased as: 
was the ratio of A^bh to A^buigo the same at all earlier 
epochs than observed now? 

The 0.3 mag displacement of the dotted line from the 
solid line in Figure 1 a implies that host galaxies at a given 
A4bh appeared brighter in the past only by (\dm/dz\) « 
0.20 mag in R. If host galaxies did only fade since z ~ 2, 
and by more than (\dm/dz\) = 0.20 mag, they would now 
appear displaced off the z = relation (solid line). If 
the host galaxies at z — 2 have bluer SED than assumed 
in our fc-corrections, the maximal fading rate would be 
even slower. For example, if we use the SED of an Scd- 
type galaxy in our fc-correction, the fading rate would be 
consistent with essentially (\dm/dz\) w mag - i.e. no 
luminosity fading. Of course, if they are to gain in A4bh, 
through accretion but not grow in A4b u ige, the allowed 
fading rate for consistency with z ~ would be smaller 
still. 

To interpret the fading rate in terms of the host stellar 
population content, we generate a grid of stellar synthesis 
models from Bruzual & Chariot (2003) with a Salpeter 
initial stellar mass function. The two sets of models 
(Fig. 2) have formation redshift zj — {2,5}, i.e. uni- 
verse age = {3.2, 1.2} Gyr, each with star formation rates 
that decrease exponentially with a decay time constant 
of t = {0,0.5,1,2,5} Gyr, where r = corresponds to 
an instantaneous burst followed by passive evolution. For 
all the models except that for r = 5 Gyr, we find that 
(\dm/dz\) « 0.9 mag over the range < z < 2. Fig- 
ure 2b illustrates how the fading scenarios translate into 
B — R color of the host galaxies. Only the models with 
t = {0, 0.5, 1, 2} Gyr evolve to early-type galaxy colors by 
z = 0; they have a fading rate (\dm/dz\) « 0.9 mag. The 
r = 5 Gyr scenario results in a present day color 0.4-0.5 
mag bluer than average early-type galaxies. 

If the z — 1...3 host galaxies were simply to fade (and 
do not gain in A4bh), resulting in present-epoch early-type 
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host galaxies we need to account for this luminosity evolu- 
tion, before comparing to the local Mbh-Mr relation. We 
do this in Figure 3 by taking the measurements of distant 
hosts (open symbols) from Figure 1 and simply shift the 
open points horizontally by an amount \dm/dz\ = 0.8 mag 
(see Figure 2), appropriate for a passively fading red bulge. 
Figures 3 a and 36 show that in this scenario, quasar hosts 
at z w 1 evolve onto the z « relationship. Some of the 
deviations may be explained by our fc-correction, or if the 
deviant galaxies underwent a starburst at z « 1. How- 
ever, for host galaxies at z « 2, what is striking is that 
nearly all would evolve to be fainter than local galaxies 
with similar mass BHs. The average offset of z « 2 RQQ 
hosts after such fading would now be at least 1.3 mag, and 
as high as ~ 2 mag, below the local Mbh-Mr relation. 
Hence our data are inconsistent with simple evolutionary 
or passively fading models for an early-type galaxy. In 
fact, the only model that is marginally consistent with the 
present-day luminosity constraint, \dm/dz\ < 0.20 mag, 
requires an evolution that has a slow decline in star for- 
mation rate of t = 5 Gyr, and that evolves to a blue 
galaxy (B — R = 1.1 mag) by z — (Fig. 26). This is only 
possible if all descendants of z ~ 2 RQQ hosts were blue 
bulges now. Even the most extreme case of an elliptical, 
fully formed by Zf = oo and passively fading, would have 
\dm/dz\ « 0.6 (van Dokkum & Franx 2001) - still three 
times larger than allowed for consistency with no mass evo- 
lution. On the other hand, for largely passive evolution of 
the host galaxy populations (with Zf < 5), the only way to 
reconcile Figure 3 with the local Mbh-Mr relationship is 
to require the hosts to undergo mergers that increase the 
galaxy stellar mass by w 3x (1.3 magnitude). However, 
in this process, they must not gain an equal proportion 
of Mbh at the merger, which only takes them parallel to 
the solid line. Our conclusions that z ~ 2 hosts appear 
too faint to be reconciled easily with present-day counter- 
parts of the same A^bh would be even stronger for any 
other plausible fc-correction one could make. 

In the discussion above we have not considered that 
Mbh may also evolve. We collect here a few other in- 
teresting implications gleaned from Figures 1 and 3: 

• If evolution occurs on the two-dimensional Mbh~ 
Mr plane, the curious fact that z = 2 hosts already 
lie close to the z = relation implies that evolution 
happens mostly along the relation, allowing for a 
modest amount of fading. Black hole growth, fading 
and the acquisition of new stars through mergers or 
star-formation seem to be balanced. 

• The BH masses for most of the host galaxies in our 
study at z > 1 are already higher than that ex- 
pected for L* galaxies (M R = —20.88, Brown ct 
al. 2001). Therefore, regardless of the stellar pop- 
ulation of the z — 2 host galaxies brighter than 
Mr< —23 mag (as in the Kukula et al. 2001 and 
Ridgway et al. 2001 sample), they cannot evolve to 
< L* galaxies today and yet agree with the local 
Mbh-Mr relation. 

• While the hosts at z > 2 may subsequently increase 
in luminosity moderately via merging, and to in- 
crease A^bh, Figure 1 shows that the most massive 



BHs may not grow by 3 to 6 times without be- 
coming more massive than the biggest BHs seen at 
low z. If the bulge mass grows proportionally with 
A4bh, this would also suggest that the host galaxies 
do not increase in mass more than a factor of 3 to 
6. 

• Figure 16 shows that the host galaxies of z = 1 
quasars arc roughly as luminous as the hosts at 
z > 2 for objects that have comparable A^bh! Fig- 
ure 1 a is somewhat equivocal on this point because 
of the small number of points, despite the smaller 
scatter than in Figure 1 6. The correspondence may 
be fortuitous. Our evolutionary scenario suggests 
that the z ~ 1 hosts are consistent with a simple 
fading scenario for an E/S0 galaxy with zj = 5, 
while the evolution of z > 2 hosts may be more 
complex. 

• Because A^BH/Albuigc ratio is higher in the past 
than today, this might suggest that the AIbh-c* 
relation has a different zeropoint earlier in time. 
Furthermore, if the dominant route of galaxy evolu- 
tion since z ~ 2 is by way of galaxy merging rather 
than in situ star formation, it would suggest that 
the AIbh-c* at high redshift has a steeper slope 
than that observed today (Gebhardt et al. 2000, 
Ferrarese & Merritt 2000). A steeper slope in the 
A4bh vs. a plot might then allow galaxies (as they 
lie in Fig. 3) to grow in mass through mergers at a 
faster rate than the black hole mass, without evolv- 
ing parallel to the local Mbh-Mr relation. 

5. DISCUSSION AND CONCLUSION 

In this study, we have tested the hypothesis that both 
the stellar bodies of the host galaxies and the central black 
holes were fully formed by Zf — 5, evolving secularly there- 
after. This is motivated by the evidence that z < 1 hosts 
of luminous quasars live in elliptical galaxies (e.g. McLeod 
& Rieke 1995, and Dunlop et al. 2003). Moreover, several 
studies (e.g. Kukula et al. 2001, Sanchez & Gonzalez- 
Serrano) argue that host galaxies as early as z > 2 may be 
fully formed, early-types, undergoing passive fading. Our 
analysis of the existing data has shown that the Mbh-Mr 
relation of high-z quasar hosts is nearly identical to that 
of low-z galaxies (Figure 1). 

We have mostly focused our attention on z « 2 objects. 
On the other hand, z « 1 hosts do appear to deviate from 
the local A4bh-ATr relationships in Figure 1, which sug- 
gests they may be fully consistent with a secularly evolving 
E/S0 model with Zf = 5. However, currently, this evidence 
is limited by small number statistics. 

The close agreement between z > 2 and z w relation 
is surprising given that the underlying correlation seen at 
z = is fundamentally one between A4bh and A^buige) 
and more superficially between AIbh and bulge luminos- 
ity. There are several ramifications of this finding, one of 
which is that the bulge luminosity does not evolve as much 
as simple stellar synthesis models would predict, thus al- 
lowing the bulge and AIbh to do so in near lock-step up- 
wards. We have shown, using a set of conservative as- 
sumptions about the stellar population of the host, that 
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the stellar bulge mass, A^buige, at a given A^bh, is prob- 
ably lower in the past than today: If there is luminosity 
evolution, but no Mbh evolution, then the Mbk/-M bulge 
ratio is larger at z > 1 by almost the same factor that the 
galaxy dims due to passive evolution. Because fading (pas- 
sive or otherwise) would cause the host galaxies at z w 2 to 
significantly over-shoot the z w relation, there has to be 
a build up in M bulge stellar mass towards z = 0, whether 
by star formation or by mergers (without increasing A^bh 
in the same proportion). This was previously suggested 
by Rix et al. (1999, 2001) and Ridgway et al. (2001). 

As BHs only increase in mass, Figure la puts a strong 
constraint on how much the host galaxies can fade to ar- 
rive at the local Mbh-Mr relation by z — 0, regardless of 
their stellar population content at z > 1. Luminous host 
galaxies [Mr< —23 mag) with luminous quasars cannot 
fade to become L* galaxies today, unless the mass esti- 
mators calibrated locally systematically and significantly 
overestimate A^bh at high z. 

While this study supplies more questions than answers 
regarding the nature of galaxy-BH evolution, we show that 
the Mbb\-Mr diagram provides very useful constraints on 
the evolutionary paths of galaxies. Thus, it may be worth- 
while to extend this diagnostic by establishing locally the 
empirical relationship between AIbh an d bulge luminos- 
ity in multiple filters, since galaxies of different morpho- 
logical types are expected to trace distinct paths in this 
space. In this parameter space of AIbh vs. multi-color 
bulge luminosity diagram the locus of points for high-z 
galaxies can be used to predict the directions in which 
galaxies evolve. Galaxies in this diagram are not allowed 
to arbitrarily roam in color and luminosity because their 
evolutionary paths must ultimately take them to the local 
A^BH-bulge luminosity relationship. To understand how 
exactly galaxies traverse the phase diagram, it would be 
useful to obtain deep imaging of the environments of high- 
z host galaxies in order to constrain their merger rates. 

Thus far, we have not discussed how the presence of dust 
may affect our conclusions. Sub-mm observations of z > 2 
quasars (Isaak et al. 2002; Bcrtoldi et al. 2003; Knudsen, 
van der Wcrf, & Jaffc 2003) indicate that at least some 
host galaxies appear to be highly star forming and hence 
presumably highly obscured. If this were generically true 
of all high-z hosts, what we observe in the rest-frame V- 
band might be the dim light that filtered through a dense 
screen. In later phases, star formation may use up gas 
and dust, while supernovae explosions may blow holes in 
the ISM, causing host galaxies to become more transpar- 
ent with time. If this were to happen, the net effect of the 
fading populations and the decreasing optical depth on the 
absolute luminosity of the host galaxies is not clear. It is 
conceivable that the host does not evolve much in opti- 
cal luminosity, if the galaxy dust extinction dropped at 
the same rate as fading of the stars. The monkey-wrench 
thrown into the interpretation by the presumption of dust 
is yet another reason why it would be worthwhile to es- 
tablish a multi-color vs. Mbh "phase" diagram, especially 
towards the rest-frame infrared. However, despite the pos- 
sibility of dust, its presence must somehow conspire with 
star formation to preserve the good correlation seen at 
high-z (Figure 1), over 4 magnitudes in host luminosity. 
This suggests that either the distribution of dust is reg- 



ular across galaxy-types, or that dust obscuration is not 
too significant in a number of galaxies by z « 2. 

In computing the fc-correction, we assumed an SED ap- 
propriate for a current-day E/S0 galaxy, but galaxies at 
z w 2 are most likely to be bluer. Therefore the hosts 
are probably fainter in the rest-frame i?-band than what 
we computed, and the Mbh-Mr relationship for z = 2 
hosts would be even more indistinguishable from z w 
normal galaxies. This conclusion relies heavily on the as- 
sumption that the techniques used to estimate A^bh at 
low z do not systematically and significantly overestimate 
A4bh at high z. Currently, all Mbh measurements are 
tied to locally calibrated values. It remains unproven that 
the AGN broad-line regions at high z have similar struc- 
ture at low z. Because understanding the link between 
AIbh and galaxy bulges is fundamentally important to a 
coherent knowledge of galaxy evolution in general, it is 
important to more fully understand how locally calibrated 
A(bh measurement techniques apply to high z. 

Lastly, we prelude in the passing that a significantly 
larger sample of quasar host galaxies at z > 2 from gravita- 
tional lensing has been analyzed to study the Mbh-Mr re- 
lationship. The results from that study further strengthen 
the conclusions presented here (Peng et al. in prepara- 
tion). 
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Table 1 

Quasar and Host Galaxy Data Compiled from Literature 



Object 


z 


Filter 




NICMOS Mag 




Quasar 
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References/Comments 












Host 
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Quasar 


Err 


Filter 


Mag 


Loud? 








(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 
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0.955 
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20.09 


0.4 


19.45 


0.3 


V* 


19.7 


N 


1, assumed quasar (B — 


V) 


=0.3 
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F110M 


20.10 


0.4 


17.90 


0.3 


V* 


19.2 


N 


1 




BVF247 
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F110M 


18.87 


0.4 


20.13 


0.3 


V* 


19.5 


N 


1 






BVF262 
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F110M 


19.85 


0.4 


19.24 


0.3 


V* 


19.4 


N 


1 






PKS 0440-00 


0.844 


F110M 


18.79 


0.4 


18.42 


0.3 


V* 


19.1 


Y 


1 






PKS 0938+18 


0.943 


F110M 


19.46 


0.4 


19.78 


0.3 


V 


18.9 


Y 


1, assumed quasar (B — 


V) 


=0.3 


3C 422 


0.942 


F110M 


18.24 


0.4 


17.85 


0.3 


V 


18.9 


Y 


1, assumed quasar (B — 


V) 


=0.3 


MRC 2112+172 


0.878 


F110M 


18.06 


0.4 


18.85 


0.3 


V 


17.8 


Y 


1, assumed quasar (B — 


V) 


=0.3 


4C 02.54 


0.976 


F110M 


19.28 


0.4 


17.57 


0.3 


V 


18.5 


Y 


1, assumed quasar (B — 


V) 


=0.3 


SGP2:11 


1.976 


F165M 


20.64 


0.75 


18.96 


0.3 


B* 


20.9 


N 


1 






SGP2:25 


1.868 


F165M 


19.88 


0.75 


19.59 


0.3 


B* 


20.7 


N 


1 






SGP2:36 


1.756 


F165M 


19.73 


0.75 


19.97 


0.3 


B* 


20.7 


N 


1 






SGP3:39 


1.964 


F165M 


19.75 
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19.53 


0.3 


B* 


20.8 


N 
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SGP4:39 
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21.59 


0.75 
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0.3 


B* 
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N 


1 






PKS 1524-13 


1.687 


F165M 


19.29 


0.75 


18.03 


0.3 


B 


20.0 


Y 


1 






B2-2156+29 


1.753 


F165M 


17.81 


0.75 


17.91 


0.3 


B 


19.7 


Y 


1 






PKS 2204-20 


1.923 


F165M 


20.63 


0.75 


18.54 


0.3 


B 


20.1 


Y 


1 






4C 45.51 


1.992 


F165M 


17.79 


0.75 


17.41 


0.3 


B 


20.1 


Y 


1 






MZZ 9744 


2.735 


F160W 


21.73 


0.5 


20.02 


0.05 


B" 


21.4 


N 


2 






MZZ 9592 
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F160W 


20.70 
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0.03 


B* 


21.8 


N 
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MZZ 1558 
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F165M 
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19.08 


0.04 


B* 


21.5 
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2 






MZZ 11408 


1.735 


F165M 
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0.4 


21.08 


0.06 


B* 


21.9 
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MZZ 4935 


1.876 


F165M 


22.00 


0.4 


21.23 


0.06 


B* 


21.8 


N 
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CTQ 414 


1.29 


F160W 


19.67 


0.1 


18.09 


0.2 


V 


19.85 


? 


3 



Note. — Col. (1): Object name. Col. (2): Rcdshift. Col. (3): HST Filter. Col. (4-7): Apparent magnitude and their published 
uncertainties, in the Vega magnitude system, corrected for extinction from Schlegcl et al. (1998). Col. (8/9): Quasar magnitude 
(corrected for extinction), in the Vega magnitude system, corresponding to the filter in Col. (8). Filters with superscript * are 
photographic magnitudes. The B and V-band magnitudes of the quasars in Kukula ct al. (2001) sample arc from Vcron-Cetty & Veron 
(1996), and references therein, while the B-band magnitudes for MZZ objects are from Zitelli et al. (1992). Where V-band magnitude is 
needed and unavailable, we used (B — V) — 0.3, which corresponds to f w oc v~ 5 . Col. (10): Radio-loud quasar or radio-quiet quasar. 
Col. (11): The photometry for each set of objects comes from the references shown. References. — (1) Kukula ct al. 2001; (2) Ridgway 
ct al. 2001; (3) Peng et al. 2005. 
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Quasar and Host Galaxy Derived Quantities 
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Fig. 1. — The relationship of the black hole mass, .Mbh, vs. bulge absolute luminosity, Mr, at low z (solid round points; compiled in 
Bettoni et al. 2003) and z > 1 (open points). Solid lines: line fitted to z solid points. Open triangles: Ridgway et al. (2001). Open 
squares and Open circles: Kukula et al. (2001). Open star: Peng et al. (2005). Panel (a): .Mbh derived using the virial mass estimate. Panel 
(b): Mbh derived assuming AiBH-luminosity relation of quasars (Peterson et al. 2004). The dotted lines fitted to the z 2, radio-quiet 
AGN, points are displaced from the z = relationship by 0.3 (a) and 0.6 (6) mag. 
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Fig. 2. — Luminosity evolution for simple star formation histories. Top: We plot the i?-band luminosity evolution of stellar populations as 
a function of redshift with formation redshifts at zj = 2 (thin lines) and zj = 5 (thick lines) for an instantaneous burst and r-modcl bursts 
(Bruzual & Chariot 2003). Bottom: We plot the colors of the objects as a function of the corresponding star formation histories. The models 
arc normalized to M|j = -20.88 (Brown et al. 2001). 
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Fig. 3. — This plot is the same as Fig. 1, but the open points are shifted horizontally by assuming that the hosts evolve passively since 
z = 5 by dMft/dz = —0.8 mag. The dotted lines fitted to the z 2, radio-quiet AGN, points are displaced from the z = relationship by 1.3 
(a) and 0.8 (b) mag. 



